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We consider the Ising model in a transverse field with long-range antiferromagnetic interactions 
that decay as a power law with their distance. We study both the phase diagram and the entangle- 
ment properties as a function of the exponent of the interaction. The phase diagram can be used as 
a guide for future experiments with trapped ions. We find two gapped phases, one dominated by the 
transverse field, exhibiting quasi long range order, and one dominated by the long range interaction, 
with long range Neel ordered ground states. We determine the location of the quantum critical 
points separating those two phases. We determine their critical exponents and central-charges. In 
the phase with quasi long range order the ground states exhibit exotic corrections to the area law 
for the entanglement entropy coexisting with gapped entanglement spectra. 



Long range (LR) interactions have attracted a lot 
of attention since they could produce interesting new 
phenomena[T}J3] . Recently there have been impressive 
advances in controlling experimentally quantum systems. 
In particular it has been shown by Britton et al. [5] 
that beryllium ions can be stored in a Penning trap, 
where an accurate laser design can induce LR Ising anti- 
ferromagnetic interactions among them. This is only the 
most recent of a series of impressive experimental results 
on using trapped ions to simulate spin models [6 - 8] . Mo- 
tivated by these results we analyze the phase diagram of 
the anti-ferromagnetic LR Ising Hamiltonian in the pres- 
ence of a transverse field (LITF). The difference with the 
standard Ising model in a transverse field (ITF) is that 
the two-body part of the Hamiltonian includes interac- 
tions among arbitrary separated pairs of spins, whose 
strength decays as a power law of their distance r, r~ a 
with a > 0. 

For the LITF we i) the determine the full phase dia- 
gram of the model as a function of a, that can be used 
as a guide for future experiments with trapped ions ii) 
quantify the increase of complexity induced by the LR 
interaction for the classical simulation of the model and 
iii) characterize the phase transitions. 

Regarding both i) and iii), we identify two different 
phases. One of them, dominated by the local part of the 
Hamiltonian, is gapped and presents patterns of quasi- 
long range order (QLRO) induced by the LR part of the 
Hamiltonian. This is exotic, since normally QLRO is as- 
sociated to gapless phases. The other, dominated by the 
LR terms of the Hamiltonian, presents anti-ferromagnetic 
LR order (LRO) in the form of Neel ground states. Be- 
tween them, we observe a line of quantum phase tran- 
sitions, whose nature depends on the value of a. They 
either are in the same universality class than the ITF 
(a > 2.25), or present new universal behaviors (for 
a < 2.25). 

Concerning ii), we focus on the entanglement entropy 
content of the ground states of the LITF as a function 
of both a and the size of the system. A common belief, 



(see, however, Ref. [§]for an updated perspective) relates 
the amount of entanglement contained in a state with its 
simulability with classical computers |10| fTTj . This prac- 
tically translates into the fact that those states that obey 
the "area law" for the entanglement, can be simulated 
classically, since their entanglement scales only with the 
area of a region rather than with its volume. In particular 
all ground states of gapped short range (SR) Hamiltoni- 
ans in ID obey the "area law" [T2|, [13] . For ground state 
of LR Hamiltonians one can expect a different scenario. 
Indeed, we show that in some cases their ground state 
still obey the area law. More interestingly, in the phase 
with QLRO, we observe unusual violations to it, in a 
gapped phase, where the local part of the Hamiltonian is 
dominant. 

Our studies complement the existing one in several 
ways. On one side most of the quantum many body lit- 
erature has focused on LR dipolar interactions decaying 
with the distance as r~ 3 |14H19| . Much less has been 
done for a generic LR interaction of the type r~ a as the 
one we consider here [2"Ul - (2"2"] . For these systems, even less 
has been done with respect to the interplay between anti- 
ferromagnetism and LR interactions. This is particularly 
interesting since the anti-ferromagnetic ITF is equivalent 
through rotation of one every two spins to the ferromag- 
netic ITF, while this is not the case for the LITF. In few 
cases the effect of LR interactions has been considered on 
top of a LRO Neel state, but the interaction considered 
was non-frustrating with respect to the LRO [23, 24]- In 
the case we consider here, all the frustration comes from 
the anti-ferromagnetic nature of the LR interaction. 

As a side result, we have improved current Matrix 
Product states (MPS) techniques. Indeed we have gener- 
alized the time dependent variational principle (TDVP) 
so that it can be used with LR interactions (alternative 
approaches can be found in the literature [151 [T6l [T51 
ITTA [55J HE|). The generalization is described in detail in 
the appendix. While the choice of using an MPS ansatz 
for LR systems could be questioned (since MPS are best 
suited for ground states of local Hamiltonians [17]), our 
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results validate this choice (see also the work |28J). In- 
deed we observe that the strongest violations to the area 
law are logarithmic in the system size as for SR criti- 
cal points where ground states can still be represented 
efficiently with MPSs [57]. 

Qualitatively, however, the logarithmic corrections 
seem to coexist with a gapped entanglement spectrum, 
a very exotic feature. Indeed our data suggest that the 
ES could present bands and gaps, even if we cannot ex- 
clude the possibility of just one gap separating the first 
eigenvalue from a continuum of them. 

The model. We study a one dimensional spin chain 
with open boundary conditions (OBC). We analyze the 
ground state of the system described by the LITF Hamil- 
tonian 

H(9, a) = sm(9) £ ir^cr* o£ + cos(0) ]T a*, (1) 
i,j 1 Jl 

where i, j are two arbitrary points of the ID chain, a > 0. 
We consider the anti-ferromagnetic phase, < 9 < ^. 
The reasons for that i) it is the interesting regime for the 
experimental results in [5] , ii) we are interested in study- 
ing the interplay between LRO and frustration; iii) the 
ferromagnetic LITF has been already studied elsewhere 

The frustration effects |16] prevent us from using stan- 
dard Quantum Monte Carlo so that we turn to matrix 
product state (MPS) techniques [SS]. We use a varia- 
tional algorithm (known as TDVP [30 ) to obtain numer- 
ically the best possible MPS for the ground state of[l] In 
order to deal with the LR, the Hamiltonian is encoded in 
a matrix product operator (MPO) [29 . This requires an 
extension of the original TDVP algorithm (described in 
the appendix). Alternative techniques are also available 

[IE1E5II25]. 

In order to establish the phase diagram and to locate 
the phase transitions we study the behavior of the entan- 
glement entropy defined as 

Sl/2 = -t r Pl/2 lo S Pl/2, (2) 

where p L / 2 = tr, 1 ...j i/2 1^) (^1 an d |^) is the ground state 
of the system. The Hamiltonian [T] has a Z 2 symmetry 
generated by Q — FL a\ . For 9 > 9 C the two body terms 
of the Hamiltonian dominate. They only commute with 
G globally so that the spectrum of the reduced density 
matrix Pl/2 (neglecting spontaneous symmetry breaking 
effects) is doubly degenerate [3T|. For 9 < 9 C on the other 
hand the local part of [T] dominates. It commutes locally 
with G so that the spectrum becomes non-degenerate. 
Close to the change of degeneracy we observe a maxi- 
mum of Sl/2 that we use as the signature for the phase 
transition. 

We then analyze the entanglement spectrum (ES) on 
both sides of the transition. It is defined in terms of the 
logarithm of the reduced density matrix 

hi = log(pi), (3) 



where pi are the eigenvalues of Pl/2- For 9 > 9 C the 
ES can be fully described by using perturbation theory 
(PT). For 9 < 9 C we observe both a perturbative and 
a non-perturbative regime for the ES depending on the 
value of a. In the non-perturbative regime we observe 
in it the appearance of bands. In the same phase the 
entanglement entropy violates the area law by exhibiting 
scaling with respect to the system size. 

Once we identify the critical point we consider the fi- 
nite size scaling of the correlation functions 

(^ /2 ^ /2+L/5 > oc L- 2A *, (^/V^ 2+L / 5 ) oc L- 2A >. 

(4) 

The corresponding exponents, as a function of a present 
two different regimes. A SR regime, where the critical 
exponents are the ones of the ITF, and a LR regime, 
where the exponents vary continuously with a. 

Numerical results We have performed several TDVP 
simulations of finite chains with length L in the range 
20 < L < 150 and OBC. The interactions encoded in 
the MPOs correctly reproduce the desired power law r~ a 
in the range of distances 1 < r < L/2 [32 . For each 
simulation, we have increased the MPS bond dimension 
X up to convergence of the ground state energy to ten 
digits when passing from one value of x to the next one 
x' (typically x' = %x)- This typically happens at values 
of x < 100. 

Phase diagram. In the anti-ferromagnetic case for all 
values of a > the system shows two phases. For val- 
ues of 9 ~ 0, 9 -C 9 c (a), the ground state |f2) can be 
understood as a perturbative modification of the prod- 
uct state locally pointing along the —1 eigenvector of 
a z . In formulas, defining a z \\) = |f), cr z \\) = - 
|fi) e=0 = Yii \lY , is independent of the value of a. This 
is a gapped phase, where elementary excitations are spin- 
flips. For values of a < 1 and 9 -C 9 c (a), the ground 
state starts to encode patterns of correlations induced 
by the LR part of the Hamiltonian that suggest the exis- 
tence of a non-perturbative regime (see Fig[3]right panel). 
However, when passing from the perturbative regime to 
the non-perturbative regime, none of the observables we 
have considered shows an anomalous behavior, so that we 
conclude that there is no sharp phase transition between 
them. 

For all the values of a we have considered, at some 
6 c (a) the system undergoes a second order phase tran- 
sition to a predominantly Neel ordered state aligned in 
the x direction. At the fixed point of the Neel phase 
(at 9 = j), the ground of the system is independent 
of a, |n) e=f = + where 

0x |+) = |+) an d <y x | — ) = — |~ )j and the square brack- 
ets indicate the elementary two-sites unit cell. The first 
excited states in these phases are kinks. The gap to them 
vanishes as a approaches zero. At a = 0, indeed, the ID 
geometry of the system is completely lost and the Neel 
state melts into an exponentially degenerate ground-state 
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Figure 1. Phase diagram of the LITF from the en- 
tanglement entropy. The half chain entanglement entropy 
provides information about the phase diagram of[l]as a func- 
tion of 9 and a. The background colors represent the value 
Sl/2 f° r a system of size L — 100 and several values of a 
in the range 0.3 < a < 3 . The maximum of it, signals the 
vicinity of a phase transition. Extrapolating its position as a 
function of L, for L = 20 • • ■ 100, we locate the position of the 
transition in the thermodynamic limit 9^°. The results are 
superimposed as solid black dots connected by a dashed line. 
The anti-ferromagnetic nature of the LR interaction favors 
the z polarized phase so that the transition always occurs on 
the right of 9 = n/4, the critical value for the ITF. 

subspace made of all possible arrangements of N/2 |+) 
states and N/2 |-) . 

The value of 9 c (a) is always larger than the one of 
the ITF transition at 7r/4. This can be understood in- 
tuitively: the slower the two body interaction decays, 
(smaller a) the more the a x part of the Hamiltonian be- 
comes frustrated. As a consequence, for small values of 
the transverse field ( large values of 9) the z polarized 
state has lower energy than the highly frustrated Neel 
state so that the system transitions to z polarized phase. 
This also explains why in the (9, a) 9 C increases with de- 
creasing a. 

The phase diagram is presented in Fig. [T] where we 
plot Sl/2 °f Eq. [2] as a function of both a and 9. For 
fixed L and a, Sl/2 nas a maximum at some given 9*. 
By extrapolating the values of 9* function of L we 
determine the location of the critical point 6™ (a). These 
points are superimposed to the colored background data 
for Sl/2 at L = 100 in black in Fig. [ljand are joined by 
a dashed line as a guide to the eye. 

In the z polarized phase, we observe two very striking 
phenomena. On one side, even if the phase is gapped, 
we observe polynomially decaying correlation functions. 
Namely (a^ /2 a x /2+r ) oc r~ a while (a z L/2 ai j/2+r ) oc r~ 2a 
(similar results were also obtained in [161 HZ]). On the 
other side we observe violations to the area law for the en- 
tanglement entropy, since of the entropy increases with- 
out saturation with the size of the blocks we have consid- 
ered. There are two different regimes for the violations 



depending on the value of a. For a < 1 we observe log- 
arithmic violations to the area law so that, by using a 
tempting analogy with the case of critical systems, (33j- 
IB"5] we can define an "effective central charge" as 

Sl/2 oc - log L. (5) 

The value we determine for c/6 are reported in the up- 
per panel of Fig. [2] They are extracted by plotting 
AS* = Sl/2 ~ Sl /2 divided by V = L/L 0l with L = 20 
(we use as a reference size to eliminate the possible con- 
stant terms in the scaling). The small dispersion of the 
curves obtained from different system sizes L around a 
single curve is a confirmation of the correctness of the 
scaling form [5] Interestingly, this effective central-charge, 
in the non-perturbative regime, varies very slowly with 
9. For a > 1 we still observe a steady growth of the 
entanglement with the size of the blocks but its behavior 
is sub-logarithmic as shown in the lower panel of Fig. [2] 
Our data for a = 3 are not conclusive. They suggest 
also there the presence of sub-logarithmic corrections for 
the sizes considered but they are so slow that we cannot 
exclude that the entropy would eventually saturate for 
larger systems. We leave this as an open issue. 

The ES of Eq. [3]can be used to distinguish between the 
perturbative and the "non-perturbative" regime in the z 
polarized phase. In the perturbative regime, (for ~ 0), 
the ES shows well defined scale separation, proportional 
to different powers of the small parameter 9. The ele- 
ments of the spectrum are dominated by the leading or- 
der at which they appear in the calculation. In the ITF 
there is a single element at each order in PT, whereas 
in the LITF ES instead, multiple eigenvalues appear at 
the same order in PT. They can be identified as parallel 
straight lines by plotting the ES in a log-log plot as a 
function of 9. The slopes of them indicate to which order 
in PT the eigenvalue belongs to, as shown in Fig. [3] right 
panel for a — 2. There we appreciate both the prolif- 
eration of eigenvalues and the wide range of validity of 
PT. We also see that the ES is dominated by eigenvalues 
appearing at most at order 9 4 in PT. 

In the same range of 9, the ES for a = 0.3 looks very 
different. In Fig. [3] right panel, we do not see neither a 
clear separations of scales, nor a well-defined power-law 
behavior of the eigenvalues with respect to 9 both foot- 
prints of the "non-perturbative" regime |36j . The eigen- 
values tend to cluster in bands (and the respective gaps) 
that are robust to changes in size (at least for the range 
of sizes we have access to) . An unsolved issue is whether 
they would survive to the thermodynamic limit. 

In both perturbative and non-perturbative regime the 
logarithmic violations to the area law coexist with a 
gapped ES. This gap is likely to survive in the thermody- 
namic limit, so that these corrections are different from 
those of a quantum critical point, where the ES gap closes 
with the system size |37| . 
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Figure 2. Violations to the area law. The entanglement 
entropy for a bipartition increases monotonically with the sys- 
tem size in the whole z polarized phase. Upper panel. 
For a < 1 the entropy scales logarithmically with the size 
of the system (small inset). The pre-factor in Eq. [5] is ex- 
tracted by plotting c/6 = ^|fr, with AS = S L/2 - S La/2 
and L' = L/L , L = 20, L = 30 ■ ■ • 100. Lower panel. For 
a = 2 the scaling is sub-logarithmic (but still monotonically 
increasing) . 
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Figure 3. Structure of the ES defined in Eq. [3]as a function 
of log(#), for < 8 < 0.36, deep in the z polarized phase. Left 
panel a = 2, the spectrum presents well separated scales 
reproducible by the lowest order in PT. Several eigenvalues 
belong to the same order in PT ( that can be distinguished as 
parallel lines ) contrary to what happens in the ITF. Right 
panel, the "non-perturbative" regime for a = 0.3 , there is 
no clear scale separation, and no clear power law dependence 
of the eigenvalues with 9. The eigenvalues tend to cluster to 
form bands separated by gaps. 



The phase transition. For the anti-ferromagnetic in- 
teraction we find a phase transition for every value of 
a > 0. This has to be compared with the ferromagnetic 
case where there is a lower critical dimension a = 1/2 
|21| . A mean field analysis around the ITF critical point 
[31 [35J suggests that the LR interactions are relevant 
for a < 2 + 2A X R driving the system to a different crit- 
ical point than the SR case, with A^ R = 1/8 being the 
scaling dimension of the a x operator for the SR ITF. For 
a = 2.25 the LR is marginal, while for a > 2.25 it be- 
comes irrelevant and one should observe the standard SR 
ITF criticality. 

We check the above statements performing a finite- 
size scaling analysis of the correlation functions [4] 
(<J x /2 ,a x /2+L/5 ) oc L' 2 ^ R . The exponents A L x R {a) 
are presented in the upper panel of Fig. [4J 2A x R (a) is 
different from 2A SR for all values of a < 2 while between 
2 and 3 it becomes very close to expected SR value 1/4. 

By studying the scaling of m Eq- [5] we can ex- 
tract the value of the central charge of the corresponding 
CFT that, for the ITF, is c = 1/2. In the whole range 
of a considered, the coefficient we obtain is systemati- 
cally bigger than 1/2. The reason for that is not clear 
but probably resides in a mixture of effects, i) the ef- 
fects of boundaries are enhanced by the LR interaction, 

ii) the system sizes we can address are still too small 
to get rid of the irrelevant contributions to the leading 
scaling [3D] (indeed our data agree with a pure logarith- 
mic scaling only for the biggest lattices L = 70 • • • 100) , 

iii) the LR could induce some marginal operator induc- 
ing corrections to the scaling difficult to control|40j . The 
corresponding plot is presented in the lower panel of Fig. 

Finally we have checked the leading power-law scaling 
of the a z correlation, an operator that already for the ITF 
is not a scaling field on its own. There we expect that 
its leading scaling is dictated by the thermal exponent 
A Z R = 1. The results for the LITF are presented in the 
central panel of Fig. [4] In the SR regime, for a > 2.25, 
the exponent we extract from the fit gives an estimate of 
the thermal exponent off by around 10% clear symptom 
of contamination with sub-leading corrections. 

Conclusions and Outlook. In this paper we have con- 
sidered the effects of a LR anti-ferromagnetic interac- 
tion on the phase diagram of the ITF, in order to both 
provide a guide to future trapped ions experiments and 
study the increase of complexity induced by the LR in- 
teractions. The resulting phase diagram shows that the 
frustration favors the z polarized phase over the x aligned 
Neel phase. For all values of a > considered we have 
located the phase transition. There we have confirmed 
that the LR interaction is relevant for a < 2.25, inducing 
critical exponents different from the ones of the ITF. We 
have determined them for the cr x a x and cr z a z correlations 
(the equivalent of the magnetic and thermal exponents 
in the SR case) . They vary continuously as a function of 
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Figure 4. Long range universality class. From a mean 
field analysis for a > 2.25 the LR become irrelevant. Upper 
panel, 2A^ R as a function of a. As expected close to a — 2 
the exponent tends to its SR value 1/4 . Middle panel, 
2Aj R does not coincide with the expected thermal exponent 
2 due to subleading corrections. Lower panel. The central 
charge extracted from Sl/2, unexpectedly, is systematically 
larger than 1/2. 



a in the range < a < 2.25. The scaling of the entan- 
glement entropy in the SR regime is used to provide an 
estimate for the central charge c of the underlying CFT, 
that turns out to be systematically larger than the ex- 
pected value 1/2. We miss a complete understanding of 
this result (that however could be a manifestation of the 
fact that the system we can address are still too small to 
see the expected asymptotic scaling) and further studies 
should be devoted to clarify it. 

It is worth to mention that the complexity of the 
ground state induced by the LR part of the Hamilto- 
nian is not significantly higher than the one of their SR 
equivalent. However we have encountered surprising vi- 
olations to the area law for the entanglement entropy, 
whose strength depends on a. The strongest violations 
are found for a < 1 and are logarithmic in the sys- 
tem size. These violations that only appear in the SR 
dominated z polarized gapped phase, seem to be always 
accompanied by a finite entanglement gap and in some 
cases by the presence of bands in the ES. Further studies 
should be devoted to check the persistence of the correc- 
tions for dipolar interactions in the thermodynamic limit. 
In the x aligned Neel phase, dominated by the LR part of 
the Hamiltonian, there are no violations to the area law. 
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Time dependent variational principle with long 
range interactions 

In oder to extract the best possible MPS description 
of the ground state of a given Hamiltonian there are sev- 
eral algorithms. Each of them present some advantages, 
i.e. the original tebd is very used given is simplicity fTU] 
while variational methods based on the energy minimiza- 
tion are often preferred since they are faster [55] . Most of 
the MPS based algorithms (including the original DMRG 
proposals |U|) rely on the fact that the MPS bond di- 
mension x grows during the computation (typically from 
X to dx where d is the dimension of the local Hilbcrt 
space), and is then reduced again to x by keeping only 
the biggest singular values of a specific bi-partition of the 
system [TU] . Recently new strategies have been developed 
based on the geometric notion of the MPS tangent plane 
[5U] that allow to optimize the MPS by solving a differ- 
ential equation, without the need of extending its bond 
dimension. Here we describe how to implement this strat- 
egy for finite chains (see also and for Hamiltonians 
encoded in MPOs . We define a state ^{A}) generated 



as a MPS from a set Aj, i = 1 ■ • • L of rank three tensors 
in the standard way. The first important part of the algo- 
rithm consists in choosing a gauge for the MPS matrices. 
We work in the isometric gauge defined in Fig. |5]i). In 
this way it is possible to give to the MPS an RG inter- 
pretation. The tensor A n at site n basically coarse grain 
the block on the right of n with the site n and project 
it to a subspace of the Hilbert space. It has dimensions 
Xm Xn+i, d and projects the tensor product Hilbert space 
built from C Xn+1 <8> d into the Hilbert space C Xn relevant 
for the description of the state. At each site of the chain 
one can define Xn+i * d — Xn MPS tangent vectors, see 
|30| . The requirement that those vectors are orthogonal 
to the original MPS vector is imposed by defining them 
through the projection onto the part of the C Xn+1 ® d dis- 
carded for the description of the original state, that has 
indeed the correct dimension Xn+i *d — Xn- From a prac- 
tical point of view, we can think of the tangent vector as 
a linear superposition of L MPS states where for each of 
them one of the original tensor A n has been replaced by 
a new tensor B n as sketched in FigjSJii) . In order to en- 
sure the orthogonality, the B n tensor are defined as the 
contraction of auxiliary tensors, (for normalization con- 
venience ) the inverse square root of the reduced density 
matrix, times a matrix of free coefficients of dimension 
X n iXn+i * d — Xn called X n , and a fixed projector V n , 
that is the ultimate responsible of the orthogonality (see 
Figj5]ii) c) ). In order to deal with the LR interactions, 
we encode the Hamiltonian in an MPO. Unfortunately 
MPOs cannot encode exactly polynomially decaying in- 
teractions, so that one needs to approximate the desired 
power law with a series of exponentials (for details see 
[25J |5B1 S3 |U). The graphical representation of the 
MPO encoding the Hamiltonian is given in Fig. |5]iii). 
If we want to obtain the ground state of a given Hamil- 
tonian, we can now start with a random MPS and solve 
the Schroedinger equation in imaginary time for very long 
times. In formula we would like to solve for \ip({A(t)})) 
the long time 

B t M{A(t)})) = -HM{Mt)}))' (6) 

A possible way to do it is to project the equation onto 
the tangent plane defined as the collection of Xn+i * d — 
Xni n = 1 ' ' ' L vectors, [45 . In formula we would like to 
find the tangent vector \T) that minimizes the distance 
from#|V>({A(i)})), 

\T*),:min m \\\T)-HM{A(t)}))\\*. (7) 

The optimal \T*) is built from a collection of {T*,n — 
1 • ■ ■ L} that are used to update the A n , 

A n (t + dt)=A n (t)-dtT*. (8) 

At the end of each step the MPS state should be 
brought back to the original isometric gauge, and the 
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procedure is iterated up to convergence. From the com- simplifications. In particular it implies that one can build 
putational point of view the definition of the tangent vec- the X n matrices directly form the Hamiltonian and the 
tors as being orthogonal to the MPS state involves several A tensors as written explicitly in Fig . [5] iv) . 
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Figure 5. The TDVP for finite chains with MPO. Geo- 
metric objects represent tensors, and lines attached to them 
represent their indexes that when connected are meant to be 
contracted. Lines interrupted by capital Greek letters repre- 
sent diagonal matrices, i) A set of rank three tensor A"'^' 1 
where n = 1 • ■ • L, a = l---\n, ft = 1, ■ • ■ x«+i> d (a), are 
required to obey the gauge fixing conditions, so that they are 
isometric (b), and in the Schmidt basis, where the reduced 
density matrix A is diagonal (c). In this way they are used 
to construct a state |^({j4})) that has a well defined coarse- 
graining interpretation (d). ii) A generic tangent vector 
\T) (a) is constructed from elementary rank three tensors B n , 
n = 1 • ■ • L (b). The requirement of its orthogonality to the 
original state (ip(A)\T) — can be obtained by building all 
Bn from a projector V n , (triangle in the figure), a matrix of 
free coefficients X n and (for convenience) the inverse of the 
square root of the reduced density matrix A n . The projector 
(c) is the one that ensures the orthogonality, by defining the 
B n in the orthogonal complement spanned by the respective 
A. iii) The Hamiltonian can be encoded in a MPO following 
the recipes in the literature [25l [26l [43l g4] . iv) The TDVP 
amounts to solving the Schrodinger equation [6] projected on 
to the tangent plane, and this implies that the X n defining 
the B n variations used to update the A n at a specific step of 
the optimization are defined as in the drawing. 



